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Abstract 

An effective interaction approach is used to describe the interactions between the spin or spin 
1/2 dark matter particle and the degrees of freedom of the standard model. This approach is 
applicable to those models in which the dark matter particles do not experience the standard- 
model interactions, e.g., hidden-sector models. We explore the effects of these effective interaction 
operators on (i) dark matter relic density, (ii) spin-independent and spin-dependent dark matter- 
nucleon scattering cross sections, (iii) cosmic antiproton and gamma ray fluxes from the galactic 
halo due to dark matter annihilation, and (iv) monojet and monophoton production plus missing 
energy at the Tevatron and the Large Hadron Collider (LHC). We combine the experimental data of 
relic density from WMAP7, spin-independent cross section from XENON100, spin-dependent cross 
section from XENON10, ZEPLIN-III, and SIMPLE, cosmic antiproton flux from PAMELA, cosmic 
gamma-ray flux from Fermi-LAT, and the monojet and monophoton data from the Tevatron and 
the LHC, to put the most comprehensive limits on each effective operator. 



I. INTRODUCTION 



The presence of cold dark matter (CDM) in our Universe is now well established by a 
number of observational experiments, especially the very precise measurement of the cosmic 
microwave background radiation in the Wilkinson Microwave Anisotropy Probe (WMAP) 
experiment The measured value of the CDM relic density is 

n CDM h 2 = 0.1126 ±0.0036 , (1) 

where h is the Hubble constant in units of 100 km/Mpc/s. Though the gravitation nature 
of the dark matter (DM) is commonly believed to be well established, its particle nature 
remains alluring except that it is nonbaryonic and to a high extent electrically neutral. 

One of the most appealing and natural CDM particle candidates is the weakly-interacting 
massive particle (WIMP). If the dark matter, generically denoted by x here, is thermally 
produced in the early Universe, the required annihilation cross section is right at the order 
of weak interaction. There may be a dynamical connection between the dark matter and 
weak-scale physics. The relation between the fractional relic density of \ relative to the 
critical density and its thermal annihilation cross section can be given by the following 
simple formula Q 

Olpb 

with (av ) being the annihilation cross section of the dark matter around the time of freeze- 
out, at which the annihilation rate could no longer catch up with the Hubble expansion rate 
of the Universe. Assuming the measured f2cDM^ 2 to be saturated by a single component 
WIMP, its annihilation cross section should be about 1 pb or 3 x 10~ 26 cm 3 s" 1 . This is 
exactly the size of the cross section that one expects from a weak interaction process, which 
implies an appreciable size of production rate of the WIMP at the Large Hadron Collider 
(LHC) as well as the event rates for direct and indirect searches that reach the sensitivities 
of dark matter experiments like XENON100 j^] and Fermi-LAT 4], || respectively. 

There have been many proposed candidates for the dark matter. Without committing 
to any particular DM model so as to perform a model independent analysis, we adopt an 
effective interaction approach to describe the interactions of the dark matter particle with 



the standard model (SM) particles. Recently, there have been a number o 
approach that deals with different observable signals in various experiments 



' works in this 



17H24J, 25M28j . One simple realization of the effective interaction approach is that the 
dark matter particle exists in a hidden sector, which communicates to the SM sector via 
a heavy degree of freedom in the connector sector. At energy scale well below this heavy 
mediator the interactions can be conveniently described by a set of effective interactions. 
The strength of each interaction depends on the nature of the dark matter particle and the 
mediator. A few models that can give rise to some of the operators in this analysis are 
described in the appendix. However, note that some popular dark matter models, such as 
supersymmetry, cannot be correctly described by this effective-interaction appraoch, because 
the dark matter particles themselves also experience the SM interactions. In this work, we 
will consider various spin nature of the dark matter particle including Dirac and Majorana 
for fermionic dark matter, as well as real and complex scalar. The most important set of 
interactions among the fermionic dark matter \ an d the light fermions / are described by 
the effective operators (xTx)(fF'f), where T and T' are general Dirac matrices contracted 
with appropriate Lorentz indices. We will discuss these and other operators in more details 
in the next section. 

One of the most anticipated signals of dark matter at hadronic colliders is a large missing 
energy in association with jets, photons, or leptons, such as monojet and monophoton plus 
large missing-energy signatures. For example, if we take one of the operators, an d 
attach a gluon or a photon to a quark leg, it will give rise to a monojet or a monophoton 
plus missing energy event. The Tevatron experiments and the LHC experiments have been 
actively searching for these signatures in some other context, such as large extra dimensions 
291 ] . We will use the most updated data on monojet and monophoton production from 



the LHC 



30| and the Tevatron 



3jJ-|33[ to constrain each effective operator. It turns out 



that the limits from the LHC and Tevatron are comparable to those obtained from indirect 

n 



detection data (PAMELA (3J, l35| and Fermi-LiAT |4|, l5j ), but inferior to those obtained from 
direct detection data (XENON100 js] and CDMS [36^) if a particular operator contributes 
;o spin-independent cross sections. Some recent works in this direction have been in Refs. 



-28 



Dark matter annihilation in the galactic halo gives rise to a number of observable signals, 
including excess in positron flux, antiproton flux, and gamma-ray over the corresponding 
cosmic backgrounds. The most current positron flux and antiproton flux data come from 



PAMELA 



34 



351 ] . The positron-fraction spectrum showed an uprising trend up to about 



100 GeV 



34 



background 



. H Nevertheless, the antiproton flux is consistent with the expected cosmic 



351 ] . The effective operators such as (xx)(qq) can gi ye r i se to dark matter 



annihilation into light quarks, which will eventually hadronize into antiprotons. In Refs. 



13 



161 ] the effects of dark matter annihilation on positron flux and antiproton flux were studied. 



It was shown that the antiproton flux data can give a better constraint on the effective dark 
matter interactions than the positron flux data 131 ] . So in this work we focus on antiproton 
flux data when we use antimatter search experiments to constrain the effective dark matter 
interactions. 

Another powerful set of indirect detection data comes from the gamma ray due to dark 
matter annihilation in the galactic halo. The data from the extragalactic sources contain 
large uncertainties, such that we concentrate on the galactic data in this work. Currently, 
the best data come from the Fermi-L AT experiment j^J. It detects gamma rays in sub- 
GeV region to hundreds of GeV from all directions (0° < |6| < 90°, 0° < I < 360°), i.e., 
including Galactic Center (GC) (0° < \b\ < 10°), low-latitude (10° < \b\ < 20°), mid-latitude 
(20° < \b\ < 60°), and high latitude (\b\ > 60°). The data on the photon spectrum from the 
low-latitude (10° < \b\ < 20°, 0° < I < 360°) [5| recorded by the Fermi-LAT indicated a 
continuous spectrum and mostly consistent with the known backgrounds. We can therefore 
use the data to constrain on additional sources of gamma-ray, namely, the annihilation of 
the dark matter into quarks, followed by fragmentation into neutral pions, which further 
decay into photons. Some recent works in using the Fermi-LAT to constrain various models 



,0 



24|. 



or effective dark matter interactions can be found in Refs. 

Another important method of detecting DM is through the direct collision between the 
DM particles in the halo with the nuclei of the detecting materials. The DM particles then 
lose a fraction of the kinetic energy to the nuclei, which can be detected by a phonon-type 
or scintillation-type or ionizat ion-type signal or some combinations of these types. Since 
the energy transfer is only of order O(10 — 100) keV and the event rate is extremely low, 
an almost background-free environment is needed. The most recent result comes from the 
XENON100 Collaboration 3|, which did not see any signal events and obtained limits on 
the spin-independent (SI) cross sections versus the DM mass. The 90% CL upper limit on 
ctsi ~ 10~ 45 cm 2 for m x = 50 GeV. We are going to use the limits presented in Ref. j^] to 
constrain the effective operators. 

1 A very recent result from Fermi-LAT 4| showed that the uprising trend continues to about 150 GeV. 
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In this work, we are going to constrain each operator from the combined data sets on 
relic density (WMAP), direct detection (XENON, ZEPLIN and SIMPLE), cosmic antipro- 
ton flux (PAMELA), cosmic gamma-ray flux (Fermi-LAT) , and monojet and monophoton 
production (Tevatron and LHC). The organization of this work is as follows. We describe 
the set of effective operators for fermionic and scalar DM and describe their nonrelativistic 
limits in Sec. II. In Sec. Ill, we calculate the relic density assuming production of DM from 
the thermal equilibrium and by solving the Boltzmann equation. In Sec. IV, we calculate 
both the spin-independent (S 
operator by the XENON100 



and spin-dependent (SD) cross sections and constrain each 



XENON10 (371, ZEPLIN (38|, and SIMPLE [39j data. In 
Sec. V and VI, we calculate the cosmic antiproton and gamma-ray flux, respectively, and 
constrain each operator from the most current data. In Sec. VII, we calculate monojet and 
monophoton production at the Tevatron and LHC, and use the most current data to con- 
strain each operator. In Sec. VIII, we perform a combined analysis by adding the chi-square 
of each data set. We conclude in Sec. IX. 

Before we close this introduction section we would spell out the improvements that are 
achieved in this work. 

1. The constraint from the LHC monojet and monophoton production was only recently 
done in Refs. 



11 



121 ] . We perform an independent analysis here. 



2. We perform a full calculation of the relic density by solving the Boltzmann equation. 



ZEPLIN 



3. In the spin- dependent cross sections, we include data sets from XENON10 [37], 
and SIMPLE 



4. We combine all data sets in the combined analysis. The resulting limits will be the 
most stringent so far. 

Based on these improvements and by demanding an operator not to give too much relic 
density to the Universe and satisfying the current experimental constraints from direct and 
indirect detection, and from collider data, a vast number of effective DM operators are 
indeed ruled out. The conclusion obtained here is important for building an effective model 
for the dark matter. 



There are recent works that consider light mediators between the SM fermions and the dark matter 



41], and light dark matter 



42 



the unitarity bound of the operators 
3 A generalized analysis of WIMP in nonrelativistic limit can be found in Ref. 



43|. 



II. EFFECTIVE DARK MATTER INTERACTIONS 



For simplicity, we will assume there is only one component of dark matter denoted by x 
and it is a standard model singlet. Here the x can stand for a Dirac or Majorana fermion, 
real or complex scalar, depending on the context. Also, / stands for a SM fermion, including 
quarks and leptons. We will include all quarks and leptons in our analysis. We briefly discuss 
in the Appendix a few hidden-sector models that can give rise to some of the operators used 
in this work in certain limits. For dark matter of spin 1 and spin 3/2, the reader may refer 



to the works in Refs. 



27 



28(]. 



The first set of operators that we consider is for fermionic DM. Its effective interactions 
with a pair of fermions include vector-, axial- vector, or tensor-type exchanges, given by the 
following dimension 6 operators 




where Aj is the heavy mass scale for the connector sector that has been integrated out and Ci 
is an effective coupling constant of order 0(1) that can be absorbed into Aj. It is understood 
that for Majorana fermion the vector and tensor structures are absent. 

6 



Next set of operators are for fermionic DM associated with (pseudo) scalar-type exchange 



07 = E ^ («)(//) , 0) 
/ 7 



Os = J2^-(x7"x)(ff) , (10) 



/ 



o 9 = J2^(xx)(h 5 f) , (11) 

/ 9 

°io = £ %^ (X7 5 X) (i¥/) • (12) 

The mj dependence in the coupling strength is included for scalar- type interactions be- 
cause this factor appears naturally from dark matter models with scalar exchange diagrams. 
Another light degree of freedom that couples to the fermionic dark matter is the gluon field 

°n = (XX) (-^-G^Gj , (13) 



Af x VAA/ V 12tt 



°* - '-w W*) (-^ G "" G -) • (") 



L 12 

C13 . /a 



O13 = (XX) (^ G " UG ^) > (15) 

°- - ^ ^ • < 16 > 

For operators involving gluons, the factor of strong coupling constant a s (2m x ) is also in- 
cluded because these operators are induced at one loop level and evaluated at the scale 2m x 
where m x is the dark matter mass. 

Finally, we also write down the corresponding operators for complex scalar DM. 

°- = E#(^)(W) , (17) 

°" = E if ( x ^ x ) (/7V/ ) ' (18) 



o 



17 



A5o 

yf 

/ 17 



E^fcW/) , (19) 



0l8 = J2 l ^(x^x)(h 5 f) , (20) 



/ 18 



o 



it VX) {-^G„) , (21) 



19 



L 19 

C20 / J- \ / OL 



°*> = t {x ' x) ^" G -y (22) 



We note that for real scalar dark matter the vector couplings in Eqs. ljTT]) and (fT8|) are absent. 
In what follows, we simply focus on the complex scalar dark matter. Note also that we have 
redefined the coefficients of some of the operators, which are different from our previous 



works 



13 



24f . such that they can conform with the normalization for the nucleon matrix 
elements used in the literature for the direct detection experiments. 

In Ref. 13[, we showed that in the calculation of the annihilation cross section for the 
DM relic density, the relative importance of each operator can be understood by considering 
the nonrelativistic expansion of the operator and studying the velocity dependence. We 
briefly review this matter here for convenience. In the nonrelativistic limit, the spinors for 
the Dirac DM \ and \ annihilation are ip ~ (£, e£) T and ip ~ (e?7^, rf)--/ where £ and rj are 
two-components Pauli spinors and e = 0(v/c). We can expand ip^^ip as 

V>7> ^ (1 + e 2 )??W 

where the spatial components are not suppressed by v/c. On the other hand, ^j^j 5 ^ in the 
nonrelativistic limit are 

^7°7 5 ^ - (1 + e 2 )V"£ 

where the spatial components are now suppressed by v/c. It is clear that in the nonrelativistic 
limit the time and spatial components of the vector and axial vector bilinear behave very 
differently. We can then consider them separately when it is contracted with the trace of the 
light fermion leg. If we look at the trace of (/7 M /) or (/7 M 7 5 /) in the annihilation amplitude, 
the time component part after being squared gives a quantity close to zero, while the spatial 
component part gives a quantity in the order of vn? x . Therefore, it is clear now that ip / ~f tJ- ip 
multiplied to (fj^f) or (/7 P 7 5 /) will not be suppressed, while -j/^^ 5 ^ multiplied to (fj^f) 
or {fl/j,1 5 f) will always be suppressed. Therefore, the operators 0\ and O3 can contribute 
to annihilation much more than the operators O2 and O4. All the other operators can be 



understood similarly 13[. From Ref. 13], we knew that some of the operators are doubly 
suppressed by the velocity of the dark matter combined with either a light fermion mass or 
strong coupling constant. Note that some of the lower limits that we obtained before are 
relative low compared to the dark matter mass. In such cases, one may question the validity 
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of the effective interaction approach. The physics behind is easy to understand. The effects 
of such operators are very suppressed because of the small velocity suppression or helicity 
suppression, not because of the size of the A. Therefore, the A has to be small enough in 
order to see an effect from these operators. We argue that the effective momentum transfer 
of such velocity-suppressed operators should be m x (v/c). With (v/c) ~ 10~ 3 for the DM 
velocity at the present epoch, as long as the ratio m x (v/c)/A remains small, we expect the 
effective interaction approach can still be valid. 

The above effective operators are relativistically invariant and therefore appropriate for 
the calculation in the relic density of the dark matter and its implication at collider physics. 
However, for direct detection experiments, we need to have a nonrelativistic reduction of 
these operators since the local dark matter velocity in the halo is of order (v/c) ~ 10 -3 . It is 
straightforward to demonstrate in the nonrelativistic limit only eight operators are relevant 
for the direct detections. These are 0\, 4 , 5 , 7 , On, Oi 5 , 0± 7 , and Oi 9 . One can further 
show that only Oi, O4 and 7 are independent, since we have the following nonrelativistic 
reduction 



5 - 


^0 4 


(23) 


On- 


^0 7 


(24) 


Ol5- 


-►Ol 


(25) 


o 17 - 


^0 7 


(26) 


o 19 - 


^0 7 


(27) 



In Table |T] we summarize some of the features of the operators discussed in this section. At 
decoupling time, v/c ~ 0.1 and hence non-relativistic reduction is no longer applicable. The 
velocity scaling behaviours for each operator shown in the last column of Table U for the 
annihilation cross sections are just merely serving the purpose to illustrate the physics. For 
our numerical work, we use the full expressions for the annihilation cross sections presented 
at the Appendix. 

In our analysis in the following sections, we will treat one operator at a time. This 
working assumption of treating one operator at one time may seem unreasonable. However 
it is a matter of choosing between controlling the number of parameters and the assumptions 
involved. If we treat each SM favor separately, then one operator at a time would mean the 
DM only couples to one quark (say u quark) but not to the other (say d quark). It would be 

9 



very strange that the new physics only couples to up quark but not to the others. But if we 
take more than one operators at the same time, the number of parameters will grow out of 
control in such an analysis. On the other hand, we have summed over all SM fermions for 
each operator. The quantum numbers of the new interaction for the SM fermions could be 
very different from one another. It is entirely model dependent. There would be too many 
parameters if we treat them all different. Even if we assumed different coefficients for each 
SM generation, we would still introduce more parameters. Here in this work, we take the 
democratic choice such that the coefficient for each SM fermion is of the same order, and we 
have treated them the same. Therefore, we sum over all SM fermions in each operator. 

We also note that the effective operators studied here in this work do not address the issue 
of gauge invariance. Imposing SU(2) gauge invariance for the SM fermions would impose 
certain relations among operators and hence their coefficients. Certain operators like those 
with an explicit factor of SM fermion mass rrif breaking SU(2) invariance explicitly can be 
made covariant by introducing the Higgs field. Such issues have been partially addressed in 
the literature, see for example in 



III. RELIC DENSITY 

In the standard cosmic picture, it is assumed that the DM particles were in thermal 
equilibrium with the other SM particles via various fundamental processes such as XX PP 
where P is any SM particles. At the high temperature Early Universe, the DM particles 
were kept in thermal equilibrium as long as the reaction rate, scaled by the temperature, 
was faster than the expansion rate H (the Hubble parameter) of the Universe. The Universe 
cooled down as it kept on expanding. At around the temperature that the reaction rate fell 
below the expansion rate H, the DM particles began to decouple from the thermal bath. 
The DM particles will keep on annihilation into the SM particles until the point that they 
could no longer effectively find one another. The remaining number density of the DM 
particles became the relic density that we can observe today. 

The input of standard computation of relic density is a ann calculated from each operator. 
For each operator, we list the (der ann / d cos 8) in the Appendix and the nonrelativistic limit 
of er anm t> in the last column of Table [I] In our numerical analysis we require the resulting 
relic density to be less than the measured value from WMAP7 [l| given by Eq. ([1]). When 

10 



TABLE I. List of properties of each operator that we define in this section. "SI" and "SD" stands 
spin-independent and spin-dependent cross sections for direct detection. 



Operator 


NR Limit 
(Direct Detection) 


SI 


SD 


Dirac / Complex 


Major ana/Real 


NR Limit (cr anni u) 
(Relic Density) 


Ox 
2 
Os 
O4 
Os 

o 6 


Yes 
No 
No 
Yes 
Yes 
No 


Yes 

No 
No 


No 

Yes 
Yes 


Yes 
Yes 
Yes 
Yes 
Yes 
Yes 


No 
Yes 
No 
Yes 
No 
No 


Nr-m 2 
ttA\ 
N c m\v 2 
67rA| 
N c m\ 
nAl 
Nctn^v 2 

6nA* 
2N c m 2 i 

nAi 
2N c m\ 
nAi 


7 
Os 

Og 
Oio 


Yes 
No 
No 
No 


Yes 
- 


No 
- 


Yes 
Yes 
Yes 
Yes 


Yes 
No 
Yes 
No 


Ncm^rn^v 2 
8?rA« 

2kA% 
Nnm 2 ,miv 2 

J A 

87rAg 
N c rn 2 m\ 
27rA? n 


On 

0x2 
Ox3 
OxA 


Yes 
No 
JNo 
No 


Yes 


No 


Yes 
Yes 
Yes 
Yes 


Yes 
No 
Yes 
No 


a 2 m^v 2 
288ir 3 A b 11 

2 4 
a s m x 

2 4 2 

1287r 3 A < i , 3 

2 4 

327r 3 A« 4 


»-> / 15 
Oi 6 


No 


1 Co 


Nn 

1 > KJ 


1 Co 

Yes 


Nn 

1> KJ 

No 


N c m 2 ,v 2 
6ttA? 5 


O17 
Oxs 


Yes 
No 


Yes 


No 


Yes 
Yes 


Yes 
Yes 


A-kA\ 7 
N c m 2 
47 rAf s 


Ox9 

o 20 


Yes 
No 


Yes 


No 


Yes 
Yes 


Yes 
Yes 


a 2 m x 
1447r 3 Af g 

4» 2 m 2 , 

30l7r 3 A| 



the relic density is smaller than the measured DM density, the DM in the Universe may 
contain more than one species. Since the relic density roughly scales inversely with <r ann , the 
WMAP7 data implies an upper limit on the heavy scale A for each operator. The results are 
shown in Fig. [1] with the requirement of the resulting density to be less than the 2a upper 
limit of the WMAP7 data. One notices that as the DM mass becomes larger all the curves 
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FIG. 1. The upper limits on A due to the WMAP7 data of Eq. ([T]). We require the resulting relic 
density less than the WMAP7 jl|: central value plus 2a error. 



increase gradually, except those for On and Oig which approach a constant. This can be 
understood by looking at the nonrelativistic limits of the annihilation cross sections listed in 
the last column of Table [H For On and 0\$, they are proportional to 1/A 4 and independent 
of m x , while for all other operators they are proportional to either m^/A 4 , m^m^/A 6 , or 
m 4 /A 6 , from which we can see that the power of A 2 in the denominator is one or two higher 
than the power of w? x in the numerator. 
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IV. DIRECT DETECTION 



The solar system moves around in the Galactic halo with a nonrelativistic velocity v ~ 
10 _3 c. When the dark matter particles move through a detector, which is usually put 
under a deep mine or a mountain to reduce backgrounds, and create collisions with the 
detector, some signals may arise in phonon-type, scintillation-type, ionization-type, or some 
combinations of them, depending on the detector materials. The event rate is extremely low 
because of the weak-interaction nature of the dark matter. There are controversies among 
various direct detection experiments. Both CoGeNT 



44| and DAMA [45| observed some 
positive signals of dark matter detection, which point to a light dark matter (~ 5 — 10 
GeV) with the crsi ~ 10~ 41 cm 2 . On the other hand, CDMS 3(| and the most recent 
XENON100 3J have found nothing and disagreed with what were found by CoGeNT and 
DAMA. In the following we will use the excluded regions of the XENQN100 data for 

!8| and SIMPLE U 



spin- independent cross sections (cr SI ), and XENON10 [37J, ZEPLIN 
data for spin-dependent cross sections (osd) versus the DM mass m x in constraining the 
effective DM interactions. 

We will be interested in the non-relativistic limit only and consider one operator at a 
time. Thus possible interference effects among different operators are ignored. 



A. Spin-Independent Cross Section 

Both Oi and O7 contribute to the spin-independent cross section. For a nuclei M with 
Z protons and (A — Z) neutrons, the cross section can be obtained as 

*&(0) = 4^IM 2 (28) 



7T 



from 0\ for Dirac DM where 



is the reduced mass and 

b N = Zb p + (A-Z)b n (30) 
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with 



&p = 2^ + T*> (31) 



A\ A 2 , 

iu rid 
^ + 2-±- 

A\ A\ 



(~iU fid 

b n = ^ + 2^. (32) 



There is no Majorana case for 0\. 
For O-j with Dirac DM, we have 



o$A0) = ^\M 2 (33) 

7T 



where 



with 



and 



f M = Zf p + (A-Z)f n (34) 

/»« = { E C ? + 1^ E ( 35 ) 

Lg=ti,(i,s Q=c,b,t J 

f<&> = 1- E /#° • ( 36 ) 

(J=M,(i,S 

For Majorana DM with the same effective operator, one should multiply the above cross 
section ( 1331) by a factor of 4. 

For On with Dirac DM, the result is the same as O7 with the following couplings 

, _ m p,n 2 f (p,n) n /„ 7 \ 

For Majorana DM, multiply the cross section by a factor of 4. 
For O15 with complex scalar, the result is 

<%r(0) = ^IM 2 (38) 

which is same as 0\ with the following replacements for the couplings in (I3T1) and (|32|) 

C M _^ C u,d ^ ^ 

Ai — ► Ai6 • (40) 
For On with complex scalar, the result is same as O7 

<&«» = ^i/.^i 2 («) 
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with fjsf = Zf p + (A — Z)f n and the following replacement in ([35 



a 



u,d 
7 

A 7 



U 17 ' 



A 



17 



For Oig with complex scalar, the result is same as O-j 



^SI (n\ f |2 

a xAfW = — A \JM\ 



An 



with fa = Zf p + (A - Z)f n and 



~~ A 3 27 G 



In our numerica 
DarkSUSY [46]. 



calculations, we will use the default values for fq P ' n ^ and f^Q 1 given in 



(42) 
(43) 

(44) 

(45) 







B. Spin-Dependent Cross Section 

For O4 with Dirac DM, its contribution to the spin-dependent cross section can be ob- 
tained as 48 1 

(46) 



«0) = ^A 2 J(J+1) 



TT 



where J is the total spin of the nuclei Af, Gf is the Fermi constant and 

A = j (ap(S'p) + a n {S n )) 



(47) 



with (iSp) and (S^) the average of the proton and neutron spins inside the nuclei respectively, 
and 



a 



E 



q=u,d,s 



1 r 19 



(48) 



with Aq^ being the fraction of the spin carried by the quark q inside the nucleon p and n. 
The following combinations of isosinglet ao and isovector a\ are often seen in the literature 



aQ — a p + a n , 
a i — a p — in ■ 



(49) 
(50) 



For a recent re-evaluation of these hadronic matrix elements using the up-to-date lattice calculation results 



of the strange quark a s term and its content in the nucleon, see Ref. 



47j]]. 
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? IG. 2. The lower limits on A due to (a) spin-independent cross section limits from XENON100 



and (b) spin-dependent cross section limits from XENON10 [37|, ZEPLIN 



and SIMPLE 



For Majorana DM with the same effective operator, one should multiply the cross section 
(126]) by a factor of 4. 

For 5 with Dirac DM, its contribution to the spin-dependent cross section is the same 
as O4 with the following replacements in (148]) 



C\ 



A, 



2Cf , 



A 



5 ■ 



(51) 
(52) 



There is no Majorana case for O5. 

The current best limits on spin-independent WIMP-nucleon cross sections come from 
XENON 100 3. InRef. Q, the collaboration searched for DM candidates in their pre-defined 
signal region, but only found 3 signal events with an expected background of 1.8 ±0.6. Based 
on that they obtained the most stringent limits on DM spin-independent elastic WIMP- 



nucleon scattering cross sections, which already exceed those of CDMS [36| in almost the 
whole mass range of m x . Therefore, we only use the XENON100 data in this analysis. 
Since the XENON100 result was presented by a 90% CL upper limit curve, we take the 
conservative choice that the central value of the WIMP-nucleon cross section for each m x 
to be zero and the la error to be the 90% CL curve divided by 1.645 (assuming a Gaussian 
distribution that 90% CL is equivalent to 1.645a.) We obtain the 2a limits on A for each 
relevant operator (note only some operators in our list can contribute to SI cross section, see 
Tabled for a summary) and show the results in Fig. 12(a). For SD WIMP-nucleon scattering 
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cross sections we use the data from XENON10 |37j, ZEPLIN [38|, and SIMPLE [39|. We 
treat the SD data in the way as how we treat the XENON100 SI data. We take the central 
value for the signal cross section to be zero and the la error for each m x is obtained by 
dividing the 90% CL curve by 1.645. We combine the chi-squares from all three experiments. 
The 2a results for A of each relevant operator that contributes to SD cross section are shown 
in Fig. EJ^b). In both SI and SD cases, we apply our formulas for the proton and neutron 
separately. Since the DM mass only enters in the SD and SI cross sections through the 
reduced mass /i X Af, which is close to the nuclei mass for large m x , one expects that 
the limits should be weaker as the DM mass grows larger following merely the constraints 
given by the experiments. This is evidently true for both SI and SD cases in Figs. [2]^a) and 
Mb), respectively. The most stringent experimental constraint for the SI case is located at 
m x w 50 GeV, while for the SD case, it is about 35 GeV. These features are also reflected 
in our figures. 

V. MONOJET AND MONOPHOTON PRODUCTION AT COLLIDERS 

In principle, dark matter particles can be directly produced in hadronic collisions. How- 
ever, it would only give rise to something missing in the detection. We therefore need 
some additional visible particles for trigger. One of the cleanest signatures is monojet or 
monophoton production, which has only a high Pt jet or p hoton balanced by a large missing 



transverse momentum. Both CDF 3J| and D0 |32|, |33[ at the Tevatron and the ATLAS 



30| at the LHC have searched for such signals, though in other context such as large extra 



dimensions. 



FIG. 3. One of the contributing Feynman diagrams for monojet or monophoton production. 

In our approach of effective DM interactions, we can attach either a gluon or a photon 
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to one of the quark legs of the relevant operators. For example, in Oi_io,i5-i8 we can attach 
a gluon or a photon line to the fermion line, such as the one shown in Fig. [31 For gluonic 
operators we can either attach a gluon line to the gluon leg or attach the whole 4-point 
diagram to a quark line such that it becomes a gg-initiated process. We then calculate 
the 2—7-3 process using FORM 49], and convolute the amplitude squared with parton 
distribution functions. The final state consists of a pair of DM particles and a gluon or a 
photon. We require the jet or photon to have a large transverse momentum according to 
the pr requirement of each experiment. 

The data sets that we used in this analysis include: (i) monojet and monophoton data 



from CDF 



3jJ, (ii) monophoton from D0 32], (iii) monojet from D0 [33], and (iv) monojet 



data from ATLAS 30]. Since the observed number of events are very close to the SM 
expectation in each experiment, we use the number of observed events and the systematic and 
statistical errors given by each experiment. For example, the observed number of monojet 
events in the very high p? selection region defined by ATLAS was 167, while the expected 
number from the SM background with the errors is 193 ± 15 ± 20 (see the last column in 



Table I of Ref. 



30].) Since the contribution from each operator does not interfere with the 



SM background, we simply add it to the SM background, and so the chi-square is 



X 



(7V(A) + N SM - N ( 



obs J 



(53) 



(15 2 + 20 2 + 167) 

where N(A) is the contribution to the event number from an operator. We show the 2a 
limits for each operator in Fig. HI We found that the chi-square is dominated by the monojet 
data of the ATLAS. Note that since the production cross section decreases as the DM mass 
increases, one expects the lower limit on the effective scale A becomes weaker. This feature 
is clearly reflected in Fig. HI 

The operators involving electrons can also give rise to monophoton events at LEP, e.g., 
(X7 M x)(e7 M e)/A 2 , by attaching an external photon line to either the electron or positron 



leg. Based on the LEP data, Refs. 



12 



50] obtained limits on the scale A. For vector-type 



interaction they obtained A > 470 — 400 GeV for m x = 10 — 80 GeV, while our limit for 
the same operator is about 800 GeV (see Fig. H(a)). Since the cross section scales as 1/A 4 
for vector-type interaction, the monophoton data from LEP would have n eglig ible effects to 



our limits. For the scalar- type interaction the limits obtained in Refs. 42|, |50J, converted to 



our convention, are about 5 GeV, while our limit is about 44 GeV (see Fig. IHb)). Again, 
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including the monophoton data from LEP would have negligible effects on our results. 



VI. INDIRECT DETECTION: GAMMA-RAY FLUX 



An important method to detect the dark matter is by measuring its annihilation products 
in Galactic halo. Current experiments can detect the positron, antiproton, gamma ray, and 
deuterium from dark matter annihilation. The Milky Way halo may contain clumps of dark 
matter, from where the annihilation of dark matter particles may give rise to large enough 
signals. There are a number of experiments dedicating to measuring the gamma-ray flux 
from DM annihilation. 

The Galactic diffuse gamma rays originate primarily from the interactions of high energy 
charged particles contained in cosmic rays with the nuclei in the interstellar medium and the 
associated radiation fields of the charged particles, via a few mechanisms briefly described 
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below. While most of them are well understood, the extra-galactic component has a larger 
uncertainty. We will choose a normalization such that the total background diffuse gamma- 
ray flux is consistent with the Fermi-LAT measurement of diffuse gamma-ray flux in the 
low-latitude. This approach is the same as the Fermi-LAT when they estimated the extra- 
galactic diffuse component {5]. 

The data on the photon spectrum from the low-latitude (10° < \b\ < 20°, 0° < / < 360°) 
j^j recorded by the Fermi-LAT indicated a continuous spectrum and mostly consistent 
with the known backgrounds. We can therefore use the data to constrain on additional 
sources of gamma-ray, namely, the annihilation of the dark matter into quarks, followed 
by fragmentation into neutral pions, which further decay into photons. The production of 
photons via neutral pions is the dominant mechanism for gamma-rays at higher energies. The 
quarks can also fragment into charged pions, which in turn decay into muons and eventually 
electrons. These electrons undergo the inverse Compton scattering and bremsstrahlung, 
which give rise to photons. The photon flux coming off inverse Compton scattering and 
bremsstrahlung tends to be dominant at lower photon energies (e.g. < 10 GeV for large 
DM mass and < 1 GeV for small DM mass). On the other hand, the synchrotron radiation 
mostly falls outside the photon energy range of the Fermi-LAT. 

The choice of the data beyond the Galactic Center is simply because the gamma-ray in 
the outside region is dominated by local sources (within our Galactic halo) and we have 
clarity in understanding the background flux and point sources within the low-latitude. On 
the other hand, the Galactic Center is supposed to have a number of known and known- 
unknown point sources, including a supermassive black hole near the Center, and perhaps 
some unknown sources too. Given the purpose of constraining the new DM interactions it 
is better to pick the data from the low-latitude region that we understand the background 
better, rather than from the Galactic Center region with less control background despite 
having a larger flux. 



A. Background Diffuse Gamma Rays 

The Galactic diffuse gamma rays originate primarily from the interactions of high energy 
charged particles contained in cosmic rays with the nuclei in the interstellar medium and the 
associated radiation fields of the charged particles, via a few of the following mechanisms. 
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(i) Gamma-rays coming from the 7r° decay, which was originated from the interactions of 
the cosmic rays with the nucleons in the interstellar medium. 

(ii) Inverse Compton scattering occurs when high energy e ± collide with the photons of 
the interstellar medium, such as CMB, star-light, and far-infrared photons. 

(iii) Bremsstrahlung photons occur when high energy e^ are deflected by the Coulomb field 
of the interstellar medium. 

(iv) Those point sources that have been identified by Fermi-LAT in the low-latitude region 



(v) Synchrotron radiation occurs when high energy e are deflected by Galactic magnetic 
field. However, synchrotron radiation only gives a very weak flux in the photon energy 
range collected by Fermi-LAT. We would not include synchrotron radiation in the 
background flux. 

(vi) An extragalactic background (EGB), which is expected to be isotropic and receives 
contributions from many sources including unresolved point sources (PS), diffuse emis- 
sion from large scale structure formation and from interactions between ultra-high 
energy cosmic ray background (CRB), and relic photons, etc. This background is 
the least determined and so a fairly large uncertainty is associated with it. Following 
Fermi-LAT we use a parameterization for the photon flux 



where A and 6 are fitted parameters (the power-law index is 7 = 5 — 2) to fit the 
extra-galactic background (EGB) by minimizing the Xbkdg fr° m the data: 



where of otal = ct^rb + a ps + ^data- Hence, we used the best-fitted point for the stan- 
dard astrophysical background (SAB). The normalization A will be varied freely when 
combining with the contribution from DM annihilation. 



15] 




(54) 




i 



$SAB(£») + $EGB(£») + gcRgjg) + $P S (£i) ~ $data(£») 




All the above sources (i) to (iii) are referred as the standard astrophysical background 
(SAB). We include the SAB, point sources, and the EGB as the background photon flux 
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in our analysis. The dominant uncertainty comes from the propa gati on parameters inside 



GALPROP 



51] . We use the best-fit model from GALPROP group 52J, in which they fitted 



to a number of isotopic ratios such as B/C, Be 10 /Be 9 , and so on. We employ the NFW 
profile with the caution that the halo uncertainty can give as much as a factor of O(10) 
change to the photon flux. Since the annihilation cross section scales as either 1/A 4 or 1/A 6 , 
this translates to the uncertainty within about 50% of the lower limit of A. 
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FIG. 5. The 2a lower limits on A for each operator due to the low-latitude gamma-ray- flux data 
ofFermi-LAT f|. 



22 



B. Dark Matter Annihilation 



In Ref. 



20[, monochromatic photon- line flux was calculated via a loop with fermions 



running in it and photons being attached to the internal fermion line. Although the photon- 
line would be a smoking-gun signal to compare with the data, the rate is suppressed because 
of the loop factor. On the other hand, photons may come from the decay of neutral pions, 
which are originated from the fragmentation of the quarks in the annihilation of the dark 
matter. The chance that an energetic quark fragments into neutral pions is high and the 
branching ratio of a neutral pion into two photons is 98.823%. Therefore, the amount of 
photons coming from the quark fragmentation is much larger than those coming off a loop 
process. Nevertheless, the spectrum of such photons is continuous and in general have no 
structure, except for a cutoff due to the mass of the dark matter. In this work, we focus 
on the continuous gamma-ray flux spectrum coming from the fragmentation of quarks into 
neutral pions, followed by their decays into photons, in the annihilation of the dark matter. 
In addition, the quarks can also fragment into charged pions, which subsequently decay 
into muons and eventually electrons. The dark matter can also directly annihilate into 
taus, muons, and electrons. The taus and muons will eventually decay into electrons. All 
these electrons undergo the inverse Compton scattering and bremsstrahlung, which may 
give rise to photons too. We include all these effects in calculating the gamma-ray flux 
from DM annihilation. Such annihilation of DM will give rise to an additional source of 
diffuse gamma-rays other than the known backgrounds. If the experimental measurement is 
consistent with the known gamma-ray background estimation, then one could use the data 
to constrain the amount of gamma-ray flux coming from the dark matter annihilation, thus 
constraining the effective interactions between the dark matter and the fermions. 

We modified DarkSUSY 46j for the effective DM interactions under consideration to 



generate the photon spectrum dN^/dE^ and e spectrum dN e ± /dE e ± of a DM annihilation 
operator of Sec. II for a particular DM mass and a selected A, say A = 300 GeV. E Note 
that the DM annihilation cross section scales as either 1/A 4 or 1/A 6 depending on operators. 
The photon and e t spectra are then fed into GALPROP with the same running parameter 

there are some PYTHIA tables, which were generated using PYTHIA 



5 Inside DarkSUSY 



53|| to 



simulate quark or gluon fragmentation into pions and kaons with a central energy (2m x ) and then de- 
cay/annihilate to 7,e + ,p, and so on. After collecting a large number of events (say ~ 10 7 ), histograms 

of dN/dE(m x , E p ) vs E p are tabulated for recycle uses. As described in the manual, the uncertainties 
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should be less than a factor of 2. 



of the best-fit model 52[. The output photon flux then includes 7r° decays, bremsstrahlung, 
and inverse-Compton scattering. This DM flux is then added to all the other astrophysical 
background fluxes and compared to the data as 



Xdm 



V-/$DM + $SAB(£ 4 ) + $EGB(£i) + $CR B (£i) + $ PS (^i) - $data(#i)\ 

? { ) (56) 



where of otal = o"crb + a ps + °data- Both the normalization A of the EGB component and 
the scale A of the DM contribution are allowed to vary freely in the fit. 

Here we adopt a simple statistical measure to quantify the effect of each DM operator. 
We calculate the 2a limit on each scale Aj while allowing the normalization A to vary, until 
we obtain a chi-square difference of Ax^ M = Xdm ~ mm (x!)M) = ^ (^ a )- For each operator 
we repeat the procedures for each DM mass. We show the results in Fig. |5j For those 
unsuppressed operators the limit is of order O(TeV). But for those operators suppressed 
by the velocity of the DM, light fermion masses or strong coupling constant, the limit is 
significantly weaker of order 0.01 — 0.1 TeV. The effects due to the onset of the heavy top 
quark in the final state are discernible by the cusps seen at some of these curves in Fig. |5j 

There is another data set on the gamma rays emitted from Dwarf spheroidal satellite 
galaxies (dSphs) of the Milky Way collected by the Fermi-L AT Collaboration {^J. They 
derived 95% CL. limits on the WIMP annihilation cross sections for a number of channels. 
The upper limits on the annihilation cross sections for the most stringent channel bb are 
(1.7 — 68) x 10~ 26 cm 3 s _1 for DM mass from 10 — 1000 GeV. If we convert these limits to 
the limits on A for the operator 0\ with the bb final state only, the limits are 0.5 — 2 TeV for 
DM mass from 10 — 1000 GeV. Including the other channels, t + t~, W + W~, would 

only lower the limits of A mildly within a factor of 2. If we compare the dSphs limits to the 
limits in our summary figures 171 to [T0l for 0±, the dSphs limits are much less stringent than 
ours. Therefore, including the dSphs data would not affect our current results significantly. 

VII. INDIRECT DETECTION: ANTIPROTON FLUX 

Annihilation of dark matter particles may give rise to large enough signals of antimatter, 
such as positron and antiproton, that can be identified by a number of antimatter search 



experiments. The most recent ones come from PAMELA 34j,|35|, which showed a spectacular 



rise in the positron spectrum but an expected spectrum for antiproton compared with their 
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?IG. 6. The 2a lower limits on A for each operator due to the antiproton flux data of PAMELA 
35]. 



backgrounds. The rise in the positron spectrum may be due to nearby pulsars or dark matter 
annihilation or decays. If it is really due to dark matter annihilation, the dark matter would 
have very strange properties, because it only gives positrons in the final products but not 
antiprotons. Here we adopt a conservative approach. We use the observed antiproton 
spectrum as a constraint on the annihilation products in XX annihilation. 

The analysis performed here is similar to that of gamma ray given in the previous section. 
We modified DarkSUSY 46( for our effective DM interactions to generate the antiproton- 
flux spectrum dNp/dTp (here conventionally the kinetic energy T is used.) The source term 
for solving the diffusion equation to obtain the antiproton spectrum is given by 

2 



v 



PCDM 

Mqdm 
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W \ dN P 



(57) 



where t] = 1/2 (1/4) for (non-)identical initial state. This source term is then fed into 



GALPROP with the running parameters of the best fit model 
is employed as in the previous gamma-ray case. 



52|. The same NFW profile 



35]. We construct the chi- 



The data set in this analysis comes from PAMELA in Ref. 
square using the 46 data points from PAMELA antiproton flux and the ratio p/p. The data 
point at the lowest energy is ignored because it does not have a central value. We included 
the solar modulation effect because it is important for the data points of the low-energy 
region. We used a modulation of 500 MV. 

We found that the background estimation of p/p using the best-fit model parameters of 
the GALPROP (which only took into account the isotopic ratios B/C, Be 10 /Be 9 , Oxygen, 
and Carbon but not the p/p ratio) fits well to the data points. We, therefore, calculate the 
2cr limit on each scale A* based on the fact that the independent background estimation 
agrees well with the data points, by a chi-square difference of A% 2 = x 2 ~~ mm (Xbk g d) = 4 
(2a). We show the resulting limits for each operator in Fig. |6j We note that the limits are 
both qualitatively and quantitatively similar to those obtained in the gamma-ray case. 

In principle, one can use the PAMELA positron spectrum to constrain the interactions. 
However, the uprising e + spectrum observed could be an indication of DM annihilation if 
there are no other known sources. In order to fit the e + spectrum the size of annihilation 
cross section a ■ v ~ 10~ 24 — 10~ 23 cm 3 s _1 , the range of which depends on the dark matter 
mass. It corresponds to Ai = 0.84 TeV for the operator 0\ when a ■ v — 5 x 10~ 24 cm 3 s" 1 
and m x = 200 GeV. Given this is the fitted value, the limit obtained would be slightly worse 
than that. It is obvious that this is negligible compared with the limit from antiproton data 
(Fig. [6]). Therefore, even including the e + data would not improve our results in a significant 
way. 



VIII. COMBINED ANALYSIS 



In this section, we do a combined chi-square analysis from all the experimental data sets 
on each effective operator. Note that the relic density from WMAP7 constrains A from 
above, while all the other experiments constrain A from below. Therefore, we combine the 
chi-squares from (i) direct detection, (ii) collider, (iii) gamma-ray, and (iv) antiproton: 

X 2 (total) = x 2 (direct) + x 2 (collider) + x 2 (gamma) + % 2 (antiproton) . (58) 
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FIG. 7. The combined analysis for 0\, O2, O3, O4, O5 and Oq. In each panel, the WMAP7 
data requires the area below the blue curve (indicated by the blue arrow) while all the other data 
requires the area above the red curve (indicated by the red arrow). The allowed region is shaded 
for 2 . 



We vary the input parameter A until the increase in chi-square is 4 units from the minimum 
value, i.e., 



A X 2 = X 2 (total) - X 2 (total) min = 4 
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FIG. 8. The combined analysis for O7, Os, Og and Oiq. In each panel, the WMAP7 data requires 
the area below the blue curve (indicated by the blue arrow) while all the other data requires the 
area above the red curve (indicated by the red arrow). The allowed region is shaded for Og. 



The limit on A thus obtained is a 2a lower limit. Together with the upper limit due to 
the WMAP7 data, we show the results for all the operators in Figs. El El El and [TOJ For 
each operator there are two curves: one from WMAP7 bounded from above and one from 
all other experimental data sets bounded from below. We indicate the allowed region by 
an arrow for each curve. Except for operators O2, Og, and 0\q, the two arrows in each 
panel are pointing away from each other, and therefore no region is allowed for all other 
operators. The working assumption here is that the effective interaction between the DM 
and SM particles thermalized the DM particles in equilibrium in the early Universe and later 
decoupled the DM particles according to the standard Boltzmann equation, and there are no 
other sources for the DM. Under this assumption most of the effective operators, except for 
O2, Og, and Oi6, cannot give a smaller interaction constrained by direct detection, indirect 
detection and collider, while at the same time provide a larger interaction allowed by the 
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FIG. 9. The combined analysis for operators On, O12, O13 and O14. In each panel, the WMAP7 
data requires the area below the blue curve (indicated by the blue arrow) while all the other 
data requires the area above the red curve (indicated by the red arrow). Allowed region for these 
operators do not exist. 



WMAP7 data. This is the main result of this work. 



Very little parameter space is allowed for most of the operators because we take the 
assumption that only one operator exists for the early universe and for present day exper- 
iments. If there are more than one operators exist at the same time, then the lower limit 
obtained by all the detection experiments (collider, indirect, and direct) will be stronger; 
on the other hand, the upper limit due to the relic density will be weaker. Therefore, there 
would be more allowable regions. Also note that if we further extends to larger masses for 
m x in almost all of the operators, there could be some allowed regions. 
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FIG. 10. The combined analysis for O15, 0\q, On, Oig, O19 and O20- I 11 each panel, the WMAP7 
data requires the area below the blue curve (indicated by the blue arrow) while all the other data 
requires the area above the red curve (indicated by the red arrow). The allowed region is shaded 
for 16 . 



IX. DISCUSSION AND CONCLUSIONS 



Since we have performed the analysis for each experimental data set and the combined 
analysis, we can easily see which data set dominates for each operator. They are summarized 
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as follows. 

• dominated by direct detection: O7, O15 

• by collider: 2 , Og, O13, O u , 16 

• by indirect detection (p and 7-ray): Oi, 3 , 4 , 5 , 6 , 8 , O10, On, 0\ 8 , O 20 

• by collider at low m x and direct detection at high m x : On 

• by collider at low m x and indirect detection at high m x : On 

• by indirect detection at low m x and direct detection at high m x : 19 

The operators 2 , Og, and 0\q that have allowed regions of parameter space are dominantly 
constrained by the collider data only. This is because these operators in the nonrelativistic 
limit (e.g. the present Universe) are highly suppressed and thus cannot contribute at any 
significant level to direct and indirect detection. 

In this work, we have exhausted all recent experimental data sets from WMAP7, di- 
rect/indirect detection, gamma-ray flux, antiproton flux and collider to obtain important 
constraints on the effective interactions of the dark matter with the SM particles. We found 
that almost all effective operators, except for 2 , Og, and 0\q, cannot give a smaller in- 
teraction that was constrained by the direct and indirect detection as well as collider data, 
while at the same time provide a large enough interaction required by the WMAP7 data 
for the relic density to avoid the Universe over closed by the DM. The result has interesting 
implications to model buildings, especially those with a heavy mediator between the dark 
sector and the SM sector. A lot of possibilities shown in this work will not work if we allow 
at the same time the new physics to give a thermal relic density and to be consistent with 
the existing data from direct and indirect detections as well as collider data. 
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Appendix A: Reductions to effective operators 



1. Darkon Model 



The darkon model 55[ consists of a real-scalar boson D, which is hidden from the SM 
interactions, except for a Higgs-portal type interaction with the SM Higgs boson 

C D = -±E D 4 - ^D 2 - \D 2 H^H . 
4 2 

After the Higgs field develops a VEV v, the interactions between the physical Higgs boson 
and the darkon are given by 

V D = ——D 4 - < + Xv2 D * _ _ \ v D 2 h , 

4 2 2 

such that the interactions of the darkon D proceed via the Higgs boson. We can write down 

the amplitude for DD — > ff as 

In the limit ^> mp, the amplitude becomes 

C = C% (ff)D 2 , 

where C/A 2 = g\v/(2mwT r n1)- This is very similar to the operator On in the case of a real 
scalar DM, with the explicit rrif dependence. 

2. Higgs Portal model for fermionic DM 

The hidden sector consists of a fermion x as the DM and a scalar boson 0, which can 
mix with the SM Higgs field H: 

£ = A 1 0xx + A 2 (0V)(# + #) . 

The amplitude for XX ~ >* J I \ after boson mixing, is given by 

£ ~ (Aixx) yr ^ 2 TT^iff) , 

[2m x y - m A h 2m w 

which becomes, in the limit m x <C m^, 

c ~ (xx) (//) , 
which is exactly the same as Oj with the explicit dependence on rrif. 
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3. Z-Z' portal model 



The hidden sector consists of a fermionic DM x an d a gauge boson Z', which then mixes 



with the SM Z boson via the kinetic mixing or Stueckelberg-type mixing [56]. The hidden 
sector does not have any SM interactions originally, but via the mixing with the SM Z 
boson, some level of interactions with the SM particles is possible. Suppose the interactions 
of the DM x an d the SM fermions are given by, after mixing and integrating out the heavy 
Z' boson, 



C = xY(9i - 9 X al")x Z, + f-figl - g f al 5 )f Z 



then the amplitude for the process XX ~~ ff can be written as 

C = Xf{£ - 9>tf)x (2mx )2- m 2 Wv - 9 f al 5 )f 

In the limit of m x <C mz, the amplitude becomes 



m z m z 

+ 9 -^4- (XY1 5 X) (/7m/) + (X7 M 7 5 X) (/7m7 5 /) 



The couplings of the Z boson to the SM fermions are all of the same order and the coupling 
to x is unknown, we can, to a crude approximation, take the overall couplings to be similar. 
Thus, we arrive at the operators 0\ to O4. Note that the coefficients for different SM 
fermions may as well be different, though they are highly model dependent, they should not 
differ from one another too much. Again, we take the crude approximation that they are 
similar and thus reduce a large number of parameters. 

A more thorough discussion on deriving various effective operators from various particle 



exchanges can be found in Ref. 



57| 
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Appendix B: Annihilation cross section formulas 



Here we list all the differential cross section formulas da^ / dz for the dark matter annihi- 
lation of the operators Oj (z =1 to 20). 

§^l^|[< + i " + 2s K-^)-8™K]. m 

d<J 3 1 N C (3f r 2 . .2 , r, / 2 2\ o 2 21 /-doN 

^7 = a| 16^ t Mm + m + (m * " m ^ " 8m - fm J ' ( 3) 

tr = ifH! [ 2( "» + *™ ) + 2s(m ' + m /> + sm w - s2 i • (B5) 

lr = Kir.?- [ 2(M ™ + ( ™' + 2s(m » + " 16m ' m ' " s2 l • (B6) 



mj N, 



c 



mj N c sfif 



G?<7 7 

(ids 

<icr 9 

dz 

dcr w 

dz ~ A% 32vr (3 X 



mj N c 
A|327r 

mj N c s/3 f 



(B7) 
(B8) 
(B9) 
(BIO) 



a 



1 



A?i 4608vr 3s2/3x ' 



a 



dan 

dz 

da\2 
dz 
dan 
dz 
dan 

~dz~ ~ Af 4 2048vr 3 J x 



A? 2 4608vr 3 (3 X 
a? 1 



A?, 2048vr3 s2/3x ' 



a 



(Bll) 

(B12) 

(B13) 

(B14) 
(B15) 
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dai 5 
dz 


1 Nq f3f . 2 . 


\ A A 

- 1) — m x + rrij 


da 16 




m f ) ) , 


dz 


da 17 


m) N c (3} 




dz 


Af 7 167i(3 x ' 




da 18 


m) N c (3 f 




dz 


Af 8 16tt/3 x ' 




da w 


a 2 1 s 





(B16) 
(B17) 
(B18) 
(B19) 
(B20) 

dz " A^o 1024vr 3 /3 X ' ^ B21 ^ 
where s, t and u are the usual Mandelstam variables, z is the cosine of scattering angle, 

u m = u — m? x — irif, t m — t — m 2 x — mj, f3 x = (1 — 4m^./s) 1 / 2 , /3f — (1 — 4mj/s) 1/ ' 2 , and 

is the color factor (3 for quarks and 1 for leptons). We have absorbed the coefficients Cj 

into Aj in these formulas. The nonrelativistic limits of OiV = Oi ■ (2/3 x ) are listed at the last 

column of Table [B 



dz Af 9 2304tt 3 p x 
da 2 o a 2 1 s 
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